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ABSTRACT 

The  problem  of  cavitating  flow  paat  a  two-dimensional  curved  obstacle  is 
considered.  Surface  tension  is  included  in  the  dynamic  boundary  condition. 
The  problem  is  solved  numerically  by  series  truncation.  Explicit  solutions 
are  presented  for  the  flow  past  a  circle.  It  is  shown  that  for  each  value  of 
the  surface  tension  different  from  zero,  there  exists  a  unique  solution  which 
leaves  the  obstacle  tangentially.  As  the  surface  tension  approaches  zero, 
this  solution  tends  to  the  classical  solution  satisfying  the  Brillouin  - 
Villat  condition.  Vanden-Broeck  considered  the  effect  of  surface  tension  on 
the  cavitating  flow  past  a  flat  plate  and  on  the  shape  of  a  jet  emerging  from 
a  reservoir.  His  results  indicate  that  the  velocity  is  infinite  at  the 
separation  points.  It  is  shown  that  these  unbounded  values  of  the  velocity 

r 

are  removed  when  the  thickness  and  finite  curvature  of  the  ends  of  the  plate 
and  of  the  ends  of  the  walls  of  the  reservoir  are  taken  into  account. 4. 
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SIGNIFICANCE  AND  EXPLANATION 


In  the  present  work  we  consider  the  effect  of  surface  tension  on  the 
cavitating  flow  past  a  curved  obstacle.  When  surface  tension  is  neglected  the 
problem  has  an  infinite  number  of  solutions.  We  show  that  this  degeneracy  is 
removed  by  solving  the  problem  with  surface  tension  and  then  taking  the  limit 
as  the  surface  tension  tends  to  zero.  Explicit  numerical  results  are 
presented  for  the  cavitating  flow  past  an  ellipse. 

The  results  presented  here  should  be  relevant  to  any  physical  situation 
in  which  the  free  surface  of  a  fluid  intersects  a  solid  surface. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


NUMERICAL  SOLUTIONS  FOR  CAVITATING  FLOW 
OF  A  FLUID  WITH  SURFACE  TENSION  PAST  A  CURVED  OBSTACLE 


Jean-Marc  Vanden-Broeck 

1.  Introduction 

We  consider  the  cavitating  flow  past  a  two-dimensional  curved  obstacle 
(see  Fig.  1).  We  neglect  the  effects  of  viscosity  and  compressibility.  The 
cavity  is  characterized  by  its  constant  pressure  pc  and  its  surface 
tension  T  while  the  fluid  has  density  p,  pressure  pc  at  infinity  and 
constant  velocity  U  at  infinity.  We  restrict  our  attention  to  obstacles 
which  are  symmetrical  with  respect  to  the  direction  of  velocity  at  infinity. 

The  problem  with  T  ■  0  has  been  considered  by  many  previous  investiga¬ 
tors  (see  Birkhoff  and  Zarantonello  [1]  for  a  review).  Efficient  numerical 
schemes  were  derived  by  Brodetsky  [2],  Vanden-Broeck  [3]  and  others.  The 
results  of  these  calculations  show  that  a  solution  exists  for  all  positions  of 
the  separation  points  A  and  B.  This  degeneracy  is  usually  removed  by 
imposing  the  Brillouin  -  Villat  condition,  which  requires  the  curvature  of  the 
free  surface  to  be  finite  at  the  separation  points  A  and  B. 

The  effect  of  surface  tension  on  the  cavitating  flow  past  a  curved 
obstacle  was  considered  by  Vanden-Broeck  [3] .  Vanden-Broeck  derived  an 
asymptotic  solution  for  T  small.  He  found  that  for  most  positions  of  the 
separation  points,  the  flow  does  not  leave  the  obstacle  tangentially.  In 
addition  he  showed  the  existence  for  each  value  of  T  of  a  particular 
position  of  the  separation  points  A  and  B  for  which  the  slope  is 
continuous  at  A  and  B.  However  he  did  not  compute  these  solutions 
explicitly. 


Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041  and 
the  National  Science  Foundation  under  Grant  No.  MCS821-5064. 


In  the  present  paper  we  solve  the  cavitating  flow  problem  numerically  by 


series  truncation.  He  assume  that  the  flow  leaves  the  obstacle  tangentially. 
Explicit  solutions  are  presented  for  the  flow  past  a  circle.  Our  results 


confirm  Vanden-Broeck's  [3]  findings.  For  each  value  of  T  there  exists  a 
unique  solution  which  leaves  the  curved  obstacle  tangentially.  As  T  ♦  0, 
these  solutions  approach  the  classical  solution  satisfying  the  Brillouin  - 
Villat  condition.  Therefore  the  degeneracy  of  the  problem  with  T  -  0  is 
removed  by  solving  the  problem  with  T  f  0  and  then  taking  the  limit  as 
T  ♦  0. 

Vanden-Broeck  [4]  investigated  the  effect  of  surface  tension  on  the 
cavitating  flow  past  a  flat  plate.  He  found  that  the  velocity  at  the 
separation  points  is  infinite  for  all  values  of  T  |  0.  He  show  that  these 
unbounded  values  of  the  velocity  are  removed  by  replacing  the  plate  by  a  thin 
ellipse. 

The  problem  is  formulated  in  Section  2.  The  numerical  scheme  is 
presented  in  Section  3  and  the  results  are  discussed  in  Section  4. 
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2.  Formulation 


Let  us  consider  the  cavitating  flow  past  a  curved  smooth  and  symmetrical 
obstacle.  We  choose  Cartesian  coordinates  such  that  the  flow  is  symmetrical 
with  respect  to  the  x-axis  (see  Fig.  1).  We  denote  by  L  a  typical  length  of 
the  obstacle  and  by  U  the  velocity  at  infinity.  We  choose  L  as  the  unit 
length  and  U  as  the  unit  velocity. 

We  introduce  the  potential  function  b$  and  the  stream  function  b|> . 

The  constant  b  is  chosen  such  that  $  =  1  at  the  separation  points  A 
and  B.  Furthermore  we  choose  $  =  0  at  x  =  y  =  0  and  *  0  on  the  free 
surface.  The  flow  configuration  in  the  complex  potential  plane  f  *  »fr+it|>  is 
illustrated  in  Fig.  2. 

We  denote  the  complex  velocity  by  ?  «  u-iv  and  we  define  the  function 
x  -  i0  by  the  relation 


u-iv 


t-ie 


vly» 


(1) 


We  shall  seek  x  -  10  as  an  analytic  of  f  *  in  the  half  plane  ♦  <  0. 

On  the  surface  of  the  cavity  the  Bernoulli  equation  and  the  pressure  jump 
due  to  surface  tension  yield  (see  Ackerberg  [5]  for  details) 


el  !i 

b  a<t> 


Jl 
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o(e 


(2) 


1  <  $  <  »  . 


Here  a  is  the  Weber  number  defined  by 


)U2L 


(3) 


We  shall  seek  solutions  for  which  the  flow  leaves  the  obstacle 


tangentially.  Therefore  we  require  the  velocity  to  be  finite  at  $  *  1. 


eT  30 

Relation  (2)  shows  that  the  curvature  of  the  free  surface  jj-  is  then  also 


finite  at  $  *  1  for  a  ^  Thus  we  impose  the  condition 

1  . 


30 

"5T  <  "  at  $ 


However  we  do  not  require  the  curvature  of  the  free  surface  to  be  equal  to  the 

curvature  of  the  obstacle.  Therefore  we  allow  a  finite  jump  discontinuity  of 
36 

3?  at  *  "  1* 

Finally  the  symmetry  of  the  problem  and  the  kinematic  condition  on  the 


obstacle  yield 


6(<J>)  =  0  t  =  0,  <j)  <  0 


F[x(<fr),y(<|>)]  -  0  ^=0,  0<4><  1 


Here  F(x,y)  =0  is  the  equation  of  the  shape  of  the  obstacle  and  8($), 
x($)  and  y($)  denote  respectively  9($,0_),  x($,0_)  and  y($,0_). 

This  completes  the  formulation  of  the  problem  of  determining  the  function 
T  -  i0.  For  each  value  of  a,  T  -  i0  must  be  analytic  in  i|»  <  0  and  satisfy 
the  boundary  conditions  ( 2 ) ,  ( 4 )  and  <  5 ) . 


3.  Numerical  procedure 

He  define  the  new  variable  t  by  the  transformation 

' -‘Sr  •  l6> 

This  transformation  maps  the  lower  half  plane  if;  <  0  onto  the  unit  circle  in 

the  complex  t-plane  so  that  the  free-surface,  the  boundary  of  the  obstacle  and 

the  negative  x-axis  go  onto  the  circumference  (see  Fig.  3). 

For  x  and  y  small  the  flow  can  be  locally  described  as  the  flow 

inside  a  right  angle  corner.  Using  (6)  we  obtain 

Vo 

C  ~  (t+1  )'2  as  t  -*■  -1  .  (7) 

At  Infinity  we  require  the  velocity  to  be  unity  in  the  x-direction.  Hence 

C  ♦  1  as  t  ♦  1  .  (8) 

We  shall  represent  C  by  the  expansion 

4  A  1/  • 

C  ■  eT~  -  (t+1)/2  l  (C  +  idn+1>tn  .  (9) 

n-0 

The  expansion  (9)  satisfies  (7).  The  unknown  coefficients  Cn  and  dR  and 

the  constant  b  have  to  be  determined  to  satisfy  (2),  (4),  (5)  and  (8). 

i  i  io 

He  use  the  notation  t  -  |t|e  so  that  the  streamline  i|>  »  0  is  given 
io 

by  t  »  e  ,  0  <  o  <  2  ir.  Using  (6)  we  rewrite  (2)  in  the  form 

5 2 -1”2 ! *T f  - 1 -  »  •<•<!  •  ,io» 

Here  t(o)  and  6(0)  denote  the  values  of  t  and  8  on  the  circumference. 

The  functions  t(o)  and  8(o)  can  be  evaluated  in  terms  of  Cn,  dn  and  b 

io 

by  substituting  t  •  e  in  ( 9 ) . 


^0  TT 

In  order  to  evaluate  —  for  0  <  a  <  — ,  we  differentiate  term  by  term 

do  2 

the  expansion  for  8(o).  This  new  expansion  can  be  expected  to  converge  to 

d0  _  it  1  r/d0i  f<J6\  1  u  . 

d?  for  0  <  o  <  -  and  to  j  w  +  l^J  ¥  ]  for  o  -  ?  where  we 

°"<a)-  °"(2)+ 


]  for  o  ■  —  where 


anticipate  a  discontinuity  in  curvature 


He  solve  the  problem  approximately  by  truncating  the  infinite  series  in 


(9)  after  N  terms •  For  convenience  we  choose  N  to  be  odd.  We  find  the 

2N  coefficients  C_  and  d_  and  the  constant  b  by  collocation.  Thus  we 

n  n 

introduce  the  2N  -  2  mesh  points 

°I  "  2N-2  (I  "  J)f  1  "  1 » • • • » 2N-2  .  (11) 

H  0 

Using  (9)  and  (11)  we  obtain  [t(o)]  ,  [0{a)l  and  [— ]  in  terms 

0=0^  0=0  ^  do  0=0  ^ 

of  the  coefficients  Cn  and  dn-  Substituting  these  expressions  into  (4)  and 

(10)  we  obtain  3(N-1)/2  equations. 

From  (1)  and  (6)  we  obtain 


dx 

do 


-b 


„  .  2  a 

2  sin  - 


— T  ~ 

e  cos  6 


(12) 


djr 

do 


-b  -T  ~ 

e  sin  6  . 


_  .20 
2  sin  - 


(13) 


dx. 


ww  dy 

Relations  (12)  and  (13)  enable  us  to  calculate  [r— 1  _  and  [— H  __  in 

do  o=o^  do  o*o ^ 


terms  of  Cn,  dn  and  b.  Using  the  trapezoidal  rule  we  then  calculate 


x{ox)  and  yfOj). 


The  boundary  condition  (5)  can  now  be  rewritten  as 


f(x(oi),  y(®j)l  "  0 


I  -  N-1 


(14) 


Relation  (14)  provides  — sj—  extra  equations.  Therefore  we  have  2N-2 

nonlinear  equations  for  the  2N+1  unknowns  Cn#  and  b.  The  condition 

(8)  provides  the  extra  equation 

[T(O)]aae0  *  0  *  d5) 

The  last  two  equations  are  obtained  by  imposing  (4)  at  o  =  and  (10)  at 

o  *  — .  Thus  we  obtain  for  each  value  of  a  a  system  of  2N+1  nonlinear 
algebraic  equations  for  the  2N+1  unknowns  Cn,  dfi  and  b.  The  system  is 


-9- 


4.  Digcueslon  of  the  results 

We  consider  the  cavitating  flow  past  a  circle.  We  choose  the  reference 
length  L  as  the  radius  of  the  circle.  Therefore  we  write 

F(x,y)  -  (x-1)2  +  y2  -  1  .  (16) 

We  used  the  numerical  scheme  of  section  3  to  compute  solutions  for 
various  values  of  a.  The  coefficients  Cn  and  dR  were  found  to  decrease 
rapidly  as  n  increases. 

Our  numerical  results  confirm  Vanden-Broeck 1 s  findings.  For  each  value 
of  a  ^  •,  there  exists  a  unique  solution  that  leaves  the  obstacle 
tangentially.  Typical  profiles  of  the  free  surface  for  a  »  0,  2  and  00  are 
shown  in  Fig.  4.  Each  of  these  profiles  is  characterized  by  a  different 
angular  position  y  of  the  separation  points.  Values  of  y  as  a  function  of 
a  1  are  shown  in  Fig.  5. 

As  a  approaches  infinity,  our  solution  tends  to  the  classical  free- 

streamline  solution  satisfying  the  Brillouin  -  Vi Hat  condition.  This 

*  0 

solution  is  characterized  by  y  »  y  ~  55  .  Therefore  the  degenerancy  of  the 

problem  with  T  -  0  is  removed  by  solving  the  problem  with  T  f  0  and  then 

taking  the  limit  as  T  -*■  0. 

0 

As  a  ♦  0,  y  ♦  90  and  the  free  surface  profile  approaches  two 
horizontal  lines  parallel  to  the  direction  of  the  velocity  at  infinity  (see 
Fig.  4). 

Similar  results  are  found  for  the  cavitating  flow  past  an  ellipse.  The 
length  L  is  chosen  as  the  half-length  of  the  axis  parallel  to  the  y-axis  and 
the  ellipses  are  parametrized  by  the  dimensionless  half-length  r  of  the  axis 
parallel  to  the  x-axis  (see  Fig.  6).  The  previous  results  for  the  flow  past  a 
circle  corresponds  to  r  *  1. 


The  limiting  configurations  a  ♦  0  and  a  ®  for  r  *  0.2  are  shown 
in  Fig.  6.  The  solution  obtained  in  the  limit  a  ♦  *  is  the  classical 
solution  satisfying  the  Brillouin  -  Villat  condition.  The  solutions  corre¬ 


sponding  to  the  limit  a  0  are  characterized  by  a  separation  point  at 
x  *  r  and  y  ■  1 . 

As  r  ♦  0,  the  ellipse  approaches  a  flat  plate  parallel  to  the  y-axis 
and  the  separation  point  of  the  limiting  solution  a  ♦  •  approaches  the 
separation  point  x  *  r,  y  ■  1  of  the  limiting  solution  a  ♦  0.  Therefore 
all  the  solutions  for  r  small  and  0  <  a  <  *  have  their  separation  points 
in  the  immediate  neighbourhood  of  x  ■  r,  y  ■  0. 

Vanden-Broeck  [4]  considered  the  effect  of  surface  tension  on  the 
cavitating  flow  past  a  flat  plate  (i.e.,  the  limiting  problem  r  *  0).  He 
provided  analytical  and  numerical  evidence  that  the  velocity  at  the  separation 
points  is  infinite  for  all  values  of  a  4  *•  Our  results  show  that  these 
unbounded  values  of  the  velocity  are  removed  by  replacing  the  flat  plate  by  a 
thin  ellipse  (i.e.,  by  solving  the  problem  with  r  4  0  and  small). 

Vanden-Broeck  [6]  investigated  the  effect  of  surface  tension  on  the  shape 
of  a  jet  emerging  from  a  reservoir.  He  also  found  unbounded  values  of  the 
velocity  at  the  separation  points.  Our  results  indicate  that  these  unbounded 
values  of  the  velocity  can  be  removed  by  taking  into  account  the  finite  thick¬ 
ness  and  finite  curvature  of  the  ends  of  the  walls  of  the  reservoir . 
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condition.  Vanden-Broeck  considered  the  effect  of  surface  tension  on  the 
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cavitating  flow  past  a  flat  plate  and  on  the  shape  of  a  jet  emerging  from  a 
reservoir.  His  results  indicate  that  the  velocity  is  infinite  at  the  separation 
points.  It  is  shown  that  these  unbounded  values  of  the  velocity  are  removed 
when  the  thickness  and  finite  curvature  of  the  ends  of  the  plate  and  of  the  ends 
of  the  walls  of  the  reservoir  are  taken  into  account. 


